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^\\ Abstract. The human sense of hearing perceives a combination of sounds 'in tune' 

if the corresponding harmonic spectra are correlated, meaning that the neuronal 
pH excitation pattern in the inner ear exhibits some kind of order. Based on this 

^ I observation it is suggested that musical instruments such as pianos can be tuned by 

minimizing the Shannon entropy of suitably preprocessed Fourier spectra. This method 

reproduces not only the correct stretch curve but also similar pitch fluctuations as in 

the case of high-quality aural tuning. 
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1. Introduction 

Western musical scales are based on the equal temperament (ET), a system of tuning 

in which adjacent notes differ by a constant frequency ratio of 2^/^^ |l|. Tuning musical 

instruments in equal temperament by ear used to be a challenging task which was carried 

out by iterating cyclically over certain interval sequences. Today this task is performed 

much more accurately with the help of electronic tuning devices which automatically 

(^ recognize the tone, measure its frequency, and display the actual pitch deviation from the 

fN theoretical value. However, if one uses such a device to tune a piano or a harpsichord 

;• exactly in equal temperament, the instrument as a whole will eventually sound as if 

it were out of tune, even though each string is tuned to the correct frequency. This 

surprising effect was first explained by O. L. Railsback in 1938, who showed that this 

perception is caused by inharmonic corrections in the overtone spectrum [2] . Professional 

aural tuners compensate this inharmonicity by small deviations, a technique known as 

stretching. The stretch depends on the specific amount of inharmonicity and can be 

visualized in a tuning chart (see Fig. [I]). 

Since the inharmonicity changes from instrument to instrument, it is difficult 
to compute the appropriate stretch by means of electronic tuning methods. Some 
appliances allow the user to preselect typical average stretches for certain classes and 
sizes of instruments. More advanced tuning devices measure the individual overtone 
spectra of all notes and compute the necessary stretch by correlating higher harmonics. 
Although the latter method gives fairly good results and is increasingly used by 
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Figure 1. Typical tuning curve of a piano (figure taken from JS]). The plot shows how 
much the fundamental frequency of each note deviates from the equal-tempered scale. 
The deviations are measured in cents, defined as 1/100 of a halftone which corresponds 
to a frequency ratio of 2^/1200 ^ 1.0005778. 



professional piano tuners, many musicians are still convinced that electronic tuning 
cannot compete with high-quality aural tuning by a skilled piano technician. This 
raises the question why aural tuning is superior to electronic methods. 

When measuring the frequencies of an aurally well-tuned piano one finds that the 
tuning curve is not smooth: rather it exhibits irregular fluctuations from note to note 
on top of the overall stretch (see Fig. fl|. At first glance one might expect that these 
fluctuations are randomly distributed and caused by the natural inaccuracy of human 
hearing. However, as we will argue in the present paper, these fluctuations are probably 
not totally random; they might instead reflect to some extent the individual irregularities 
in the overtone spectra of the respective instrument and thus could play an essential 
role in high-quality tuning. Apparently our ear can find a better compromise in the 
highly complex space of spectral lines than most electronic tuning devices can do. 

As a first step towards a better understanding of these fluctuations, we suggest that 
a musical instrument can be tuned by minimizing an appropriate entropy functional. 
This hypothesis anticipates that a complex sound is perceived as 'pleasant', 'harmonic' or 
'in tune' if the corresponding neuronal activity is ordered in such a way that the Shannon 
entropy of the excitation pattern is minimal. The hope is that such an entropy-based 
optimization allows one to find a better compromise between slightly detuned harmonics 
than a direct comparison of selected spectral lines. 
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2. Harmonic spectrum, musical scales, and temperaments 



Sound waves produced by musical instruments 
involve many Fourier components. The simp- 
lest example is the spectrum of a vibrating 
string [5]. Depending on the excitation 
mechanism, one finds not only the fundamental 
mode with the frequency /i but also a large 
number of higher partials (overtones) with 
frequencies f2, fs, . . . ,. For an ideal string 
the frequencies of the higher partials are just 
multiples of the fundamental frequency, i.e. 




fn = nfi. 



(1) 



Such a linearly organized spectrum of over- 
tones is called harmonic. 



Figure 2. 

Harmonic modes of a string [4] . 



Since harmonic overtones are ubiquitous in Nature, our sense of hearing prefers 
intervals with simple frequency ratios for which the spectra of overtones partially 
coincide. Examples are the octave (2:1), the perfect fifth (3:2), and the perfect fourth 
(4:3), which play an important role in any kind of music. On the other hand, music 
is usually based on scales of notes arranged in octave-repeating patterns. Since the 
frequency doubles from octave to octave, it grows exponentially with the index of the 
notes. This exponentially organized structure of octave-repeating notes is in immediate 
conflict with the linear spectrum of the harmonics. A musical scale can be seen as 
the attempt to reconcile these conflicting schemes, deflning the frequencies in such a 
way that the harmonics of a given note coincide as much as possible with other notes 
of the scale in higher octaves. As demonstrated in Fig. [3} this leads quite naturally to 
heptatonic scales with seven tones per octave, on which most musical cultures are based. 
In traditional Western music the seven tones (the white piano keys) are supplemented by 
flve halftones (black keys), dividing an octave in twelve approximately equal intervals. 

As the twelve intervals establish a compromise between the arithmetically ordered 
harmonics and the exponentially organized musical scale, their sizes are not uniquely 
given but may vary in some range. Over the centuries this freedom has led to 
the development of various tuning schemes, known as intonations or temperaments, 
which approximate the harmonic series to a different extent. One extreme is the 
just intonation, which is entirely built on simple rational frequency ratios. As 
shown in Fig. |3)d, the just intonation shares many spectral lines with a suitable 
harmonic spectrum. However, this tuning scheme is not equidistant on a logarithmic 
representation, breaking translational invariance under key shifts (transpositions). 
Therefore, the just intonation is in tune only with respect to a speciflc musical key 
(e.g. C Major), while it is out of tune in all other musical keys. 
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Figure 3. Harmonic overtone spectrum compared with musical temperaments 
(logarithmic scale). Lower part: (a) Fundamental frequency /i = 11 Hz and the 
corresponding series of harmonics, (b) Just intonation in C-Major with octave- 
repeating patterns of non-equidistant pitches, (c) Equally tempered intonation 
with constant pitch differences. The upper part of the figure shows a zoom of 
the octave C4-C5. As can be seen, the heptatonic scale (the white piano keys, 
thick lines) of the just intonation (b) matches perfectly with the harmonics in 
(a), while the half tones (black keys, thin lines) do not lock in. Contrarily the 
equal temperament (c) deviates in all tones except A440 but it is equidistant 
and therefore invariant under shifts (transpositions) of the musical key. 



With the increasing complexity of Western music and the development of advanced 
keyboard instruments such as harpsichords, organs and pianos, more flexible intonations 
were needed, where the musical key can be changed without renewed tuning. Searching 
for a better compromise between purity (rational frequency ratios) and temperament 
independence (transposition invariance) various tuning schemes have been developed, 
including the famous meantone temperament in the renaissance and the well-tempered 
intonation of the baroque era. Since the 19*^ century Western music is predominantly 
based on the aforementioned equal temperament, which is fully invariant under key 
shifts. In equal temperament, frequencies of neighboring tones differ by the irrational 
factor 2^'^^ so that they are equidistant in a logarithmic representation (see Fig. pp). 
However, this invariance under key shifts comes at the price that all intervals (except 
the octave) are slightly out of tune, but apparently our civilization learned to tolerate 
these discrepancies. 
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Figure 4. Left: Inharmonicity coefficients B of an upright piano. The two parts 
of the data correspond to the two diagonahy crossed bass and treble sections of 
the strings. Right: Power spectrum of the leftmost string. The red arrow marks 
the fundamental frequency of 27.5 Hz. The blue arrows indicate particularly 
weak partials which are suppressed due to the position of the hammer. 



3. Inharmonicity 



The harmonic series of overtones /„ = nfi is valid only for ideal oscillators whose 
evolution is governed by a linear second-order partial differential equation. In realistic 
musical instruments there are higher-order corrections in the force law which lead 
to small deviations from the harmonic spectrum. The degree of inharmonicity is 
characteristic of each instrument and accounts for much of the color and texture of 
its sound. 

Inharmonicity in string instruments is caused by the circumstance that a realistic 
string is an intermediate between an ideal string and a stiff bar. An ideal string vibrates 
according to the differential equation y oc —y" with a linear dispersion / oc \k\, while a 
stiff bar is known to evolve according to a fourth-order differential equation y oc —y"" 
with a quadratic dispersion / oc fc^. Therefore, the stiffness of realistic strings causes 
lowest-order corrections of the form 



y oc-y 



ey 



foce + ek' 



so that the spectrum of a string is given by 



fn oc nf 1^/1 + Brfi 



n = 1,2, ... 



(2) 



(3) 



where /i denotes the fundamental frequency and /„ is the frequency of the n^^ partial. 
The dimensionless number B is the so-called inharmonicity coefficient which depends 
on the length, diameter, tension and material properties of the string. In a piano the 
value of B varies typically between 0.0002 for bass strings up to 0.4 for treble strings 
(see left panel in Fig. |4J). Strong inharmonicity can be recognized as an unpleasant 
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Figure 5. Harmonic spectrum of partials of an octave in a linear representation. 
The octave is perceived as pleasant since every second partial of the lower tone 
locks in with one of the a partials of the upper one. 



dulcimer-like sound. An important part in the art of piano construction is to keep the 
inharmonicity as small as possible. 

4. Perception of being in tune 

As mentioned before, intervals with simple rational frequency ratios are perceived as 
pleasant. In this context it is important to note that the human ear, when hearing two 
different tones simultaneously, cannot evaluate the frequency ratio of the fundamental 
modes directly: rather it recognizes coincidences in the corresponding overtone spectra. 
For example, if we hear an octave, say A2-A3, our ear compares the partials 2,4,6,... of 
the lower tone with the partials 1,2,3... of the upper one, perceiving the octave as being 
'in tune' when the two harmonic series lock in (see Fig. p«. 

In the presence of inharmonic corrections, however, it is no longer possible to 
match the two series exactly. In this case our ears search for the best possible 
compromise, minimizing the frequency differences between almost coinciding low-lying 
partials. These tiny differences are heard as so-called beats, a superposition of enveloping 
modulations of a few Hertz which an aural tuner tries to make as slow as possible. As 
demonstrated in Fig. |6j in an octave such a compromise can be achieved by slightly 
increasing the pitch of the upper tone. This means that we perceive an interval as 
correctly tuned if it is slightly out of tune in the mathematical sense. This correction, 
called stretch, plays a major role in the practice of tuning, even when the inharmonicity 
of the instrument is small. 

Today high-priced electronic tuning systems are available which can compute the 
appropriate stretch for individual instruments. To this end some of the tones are 
recorded and the inharmonicity coefficients of the strings are estimated by identifying 
the partials in the corresponding spectra. The stretch is then computed e.g. by selecting 
sequences of octaves and stretching them in such a way that the fourth partial of the 
lower tone coincides with the second partial of the upper (4:2-tuning). More specifically, 
enumerating the piano keys hjk = l...K and denoting by fn the n}^ partial of the 
k^^ string, this 4:2 method provides K — 12 equations for the octave stretches of the 
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form 



n 



(fc+12) 



r 



(k) 



n 



(k) 



(fc+12) 



(4) 



r 



where r^ = /« //i is the ratio of the n^^ partial to the fundamental frequency of the 
string. By taking the logarithm this turns into a system of i^— 12 linear equations for the 
K unknown fundamental frequencies /} . The remaining 12 unknowns are determined 
by the reference pitch A440 and the choice of the temperament. Equal temperament can 
be approximated by introducing a quadratic penalty function for the change of adjacent 
interval ratios. By solving these equations one can translate the measured partials 
directly into a tuning curve. If the inharmonicity coefficient is a piecewise smooth 
function (like the one shown in Fig. H|, the tuning curve will be piecewise smooth as 
well. Likewise one can use a 6:3 tuning scheme, which produces an even larger stretch. 
The overall magnitude of the stretch is therefore not strictly defined but rather a matter 
of taste. Some devices even interpolate between 4:2 and 6:3 stretching in order to get a 
more acceptable compromise. 

Computing the stretch by direct comparison of partials as described above yields 
piecewise smooth tuning curves. However, as mentioned in the Introduction, aural 
tuners produce tuning curves with pronounced fluctuations on top of the overall stretch, 
especially in the bass and in the treble. One of the main messages of this work is the 
conjecture that these fluctuations are not random but to some extent essential for a 
good tuning result. 

The fluctuations may have different reasons. On the one hand, each partial couples 
differently to the resonator of the instrument (the soundboard of the piano), leading 
to additional frequency shifts so that the inharmonic spectrum deviates slightly from 
the predicted form in Eq. ([3]). Another reason is the highly irregular intensity of the 
partials. As shown in the right panel of Fig. |4| the spectrum of a piano string may 
consist of dozens of partials, but even adjacent partials may differ in their power by 
more than a magnitude. Even worse, some of the partials (indicated by blue arrows 
in the flgure) are strongly suppressed if the hammer hits the string at a node of the 
corresponding vibrational mode. This suggest that in realistic situations the perception 
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Figure 6. Compensating inharmonicity by stretching octaves (see text). 
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of being in tune does not only depend on the frequency of the partials but also on their 
amplitude. 



5. Psychoacoustic aspects 

As tuning can be understood as the search for a compromise in matching higher partials, 
it will significantly depend on the acoustic and psychoacoustic properties of the inner 
ear. Psychoacoustics is a research field on its own (see e.g. [6]-[8]) and plays an important 
role e.g. in lossy data compression methods such as MPS. Here we only sketch a few 
basic elements which are essential for the method presented below. 

Let us first consider the frequency range of the ear. Starting point is a sound wave, 
which can be described as a time-dependent pressure variation p(t). Its complex- valued 
Fourier transform is given by 

p(/) = -^ fdte'-^f'pit), (5) 

where p{—f) = p*{f)- The corresponding power spectrum 

/(/) = iPim (6) 

describes the energy density of the spectral line at frequency /. As a technically useful 
measure one defines the logarithmic sound pressure level (SPL) 

^(/) = 101ogio(^) (7) 

measured in decibels (dB), where Iq refers to the hearing threshold. 

Depending on the frequency the SPL will be correlated with a certain mechanical 
response in the inner ear. Since the physical transmission mechanism is highly complex, 
one usually approximates this relationship by certain weighting functions. Below 55 dB 
the most commonly used one is the so-called A-weighting according to the international 
standard lEC 61672:2003 with the filter function 

12200^ f'* 

RaU) = ^ (8) 

(/2 + 20.62)(/2 + 122002) ^(/2 + 107.72) {p + 737.92) 

which defines the A-weighted sound pressure level (SPLA) 

LAif) = (2.0 + 20 logio RaU)) L{f) (9) 

in units of A-weighted decibels (dBA). The SPLA can be considered as a rough measure 
of frequency-dependent energy deposition in the cochlea. 

The receptor cells in the inner ear convert the excitation pattern into a certain 
neuronal response which causes the auditory perception in the brain. The neuronal 
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Figure 7. Shannon entropy as a measure for the coincidence of spectral lines. 
The figure show the superposition of two Gaussian functions representing two 
partials. If the two partials are sufficiently separated, the (continuous) entropy 
H = — J_ f{x) log2(/(x))d2; gives a constant value H ~ 4.094. When the two 
partials begin to overlap (audible as heats) ^ the entropy decreases and reaches 
a minimum {H ~ 2.094 in this example) if they coincide. 



processing is even more complex and not yet entirely understood. For this reason one 
uses a psychoacoustic measure for the perceived intensity, the so-called loudness N{f), 
which is an empirical psychological quantity averaged over many test persons. According 
to the literature this relationship is well approximated by a piecewise exponential 
function and a power law: 



N{f) 



2(iA(/)-40)/10 

(L^(/)/40)2-86 



if LAif) > 40dBA 
if LaU) < 40dBA 



(10) 



Not only the sensitivity of the ear is frequency-dependent but also its ability to 
discriminate between different frequencies. In the literature different measures for the 
frequency resolution are reported, of which the so-called just noticeable difference (jnd) 
plays the role of a lower bound [6] . The jnd is usually approximated by 



A/ 



3 Hz 
0.006/ 



if / < 500 Hz 
if / > 500 Hz . 
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6. Entropy-based tuning scheme 



We now suggest a simple entropy-based tuning scheme for musical instruments. It 
is motivated by the observation that tuning can be understood as the search for the 
best possible compromise in matching higher partials, and the main idea is that this 
compromise is characterized by a local minimum of the entropy of the intensity spectrum. 
This is highly plausible since the entropy of two spectral lines decreases as they begin 
to overlap (see Fig. [T]). 
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To test this idea we individu- 
ally recorded all keys of an aurally 
tuned piano, computed their power 
spectra and reorganized them in 
logarithmic bins in order to account 
for the finite frequency resolution of 
the inner ear. Furthermore, we re- 
moved the pitch differences, reset- 
ting the tune of all tones to equal 
temperament. Then we applied the 
following simple zero-temperature 
Monte-Carlo scheme of statistical 
physics (see Fig. |8| technical details 
are given in the appendix): 



Choose a random key and 
change its pitch randomly 
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Figure 8. Monte-Carlo scheme. 



• Add the A-weighted power spectra of all 88 tones and compute the entropy. 

• Randomly change one the of the pitches and compute the entropy again. 

• If the entropy is lower accept the pitch change, otherwise restore the previous value. 

This simple procedure is iterated until no further improvement is obtained, meaning 
that the algorithm has found a local minimum of the entropy. Note that by adding up 
all tones, the method is inherently sensitive to all intervals, not only to octaves. 




40 
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Figure 9. Typical result of the tuning procedure described in Section 6 (red 
curve) compared with the tuning curve produced by aural tuning (black curve). 
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7. Discussion 

Fig. [9] shows the resulting tuning curve of a typical run compared with the actual curve 
produced by an aural tuner for an upright piano. As can be seen, not only the overall 
stretch is predicted correctly but even the fluctuations of the two curves are highly 
correlated, especially in the bass and the treble. Apparently the entropy-based tuning 
method is capable of generating the same individual deviations from the average stretch 
as an aural tuner. This is surprising and not yet understood, but it indicates that these 
fluctuations are reproducible and may play an essential role in the practice of tuning. 

The implementation of the method is very easy. The tones are recorded, Fourier- 
transformed, mapped pointwise by the psychoacoustic filtering functions as described 
above, binned logarithmically, added up, and finally plugged into the entropy functional. 
An explicit identification of higher partials and the measurement of the inharmonicity 
is not needed. The method is expected to take automatically any anomalous spectral 
properties of the instrument into account. 

However, the method suggested here is still in an immature state. It could be 
modified in various respects and a systematic study is still outstanding. Moreover, the 
method was tested so far with only one instrument. The main open questions are the 
following: 

• Apparently there are many local minima, so that the algorithm outlined above 
produces similar but not reproducible results. 

• The Monte-Carlo results presented above were based on the A-weighted spectra 
(SPLA) in Eq. ^. If one uses instead the loudness defined in Eq. (10), one obtains 
unreasonably stretched tuning curves in the bass. 

• The spectra were logarithmically binned in units of one cent. This models a 
frequency resolution of one cent, which is smaller than the just notable difference 
{jnd). However, convolving the spectra with a frequency-dependent Gaussian 



according to the expected jnd in Eq. (11) does not improve the results. 



• More advanced Monte Carlo techniques such as simulated annealing have not yet 
been tested. 

• Instead of adding up the spectra of all piano keys, we tried to work with subsets of 
octaves, fifths and fourths, imitating the practice of aural tuners. This destabilizes 
the method, probably driving the pitches out of equal-tempered into just intonation. 
Apparently the summation over all keys allows the system as a whole to stay in 
equal temperament. 

Regarding possible technological realizations, it could be interesting to develop electronic 
tuning devices using a hybrid method, which first compute an approximate tuning curve 
by matching higher partials, and then optimize the fiuctuations of the pitches by search- 
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ing for a suitable local minimum in the vicinity. 
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Appendix A. Technical details 

Data recording and preprocessing 
(i) Record the piano keys k = 1 . . . K in WAV format. Extract the binary PCM 

(k) 

amplitudes and convert them to a series of floating point numbers i/j gM. with an 
index j = . . . ST — 1, where S = 44100 Hz is the sample rate and T ^ 20s is the 
recording time. 

(ii) Apply a fast Fourier transform (e.g. package f f tw3) to obtain the spectra yqEC 
indexed by q = . . .Q, where Q = ST/2 (the other half of the data is complex 
conjugate). The q^^ component corresponds to the frequency f{q) = q/T. 

(iii) For each k coarse-grain the power spectrum \yq p G M"*" by logarithmic binning. 
To this end define an array Im € M"*" corresponding to the frequencies f{m) = 
10 ■ 2"^/i200 g2 with m running from zero (10 Hz) to 12000 (10 kHz). Let 

Q 

4'^ := Yl <[i200+ios,(^)] |yf f , (A.l) 

9=0 

where [■] denotes rounding to an integer. Note that in this representation adjacent 
bins differ by a frequency ratio of one cent. 

(iv) Map the intensities Im to the corresponding A- weighted sound pressure levels 
(SPLA) L^). 

(v) After this preprocessing the pitch of a tone can be increased or lowered by c cents 
through a shift of the array index m — )■ m — c. This allows us to tune the instrument 
virtually on the computer. Remove the recorded pitch deviations by tuning the first 
partial to equal temperament, i.e. /} = 440 ■ 2('^~'^o^/^^ Hz, where ko is the index 
of A440. This means that the recorded stretch corrections are initially eliminated. 

Monte- Carlo dynamics 

(i) Change one of the pitch differences randomly by ±1 cent, 
(ii) Compute the sum pm = X]fc=i -^™ '-'^ ^^^ SPLA over all keys, 
(iii) Normalize Pm such that X^^Pm = 1- 
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(iv) Compute the Shannon entropy H = — J2mPm^^Pm- 
(v) If the entropy decreases, keep the pitch change, otherwise restore the old pitch. 

This procedure is repeated until no further changes take place. 
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